Bound magnetic polarons in the very dilute regime 
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We study bound magnetic polarons (BMP) in a very diluted magnetic semiconductor Cdi_ x Mn x Te 
(a; < 0.01) by means of site selective spectroscopy. In zero magnetic field we detect a broad and 
asymmetric band with a characteristic spectral width of about 5 meV. When external magnetic fields 
are applied a new line appears in the emission spectrum. Remarkably the spectral width of this line 
is reduced greatly down to 240 peV. We attribute such unusual behavior to the formation of BMP, 
effected by sizable fluctuations of local magnetic moments. The modifications of the optical spectra 
have been simulated by the Monte-Carlo method and calculated within an approach considering the 
nearest Mn ion. A quantitative agreement with the experiment is achieved without use of fitting 
parameters. It is demonstrated that the low-energy part of the emission spectra originates from the 
energetic relaxation of a complex consisting of a hole and its nearest Mn ion. It is also shown that 
the contribution to the narrow line arises from the remote Mn ions. 

PACS numbers: 75.50.Pp, 78.20.Ls, 75.30.Hx 
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An intriguing phenomenon in diluted magnetic semi- 
conductors (DMS) is the formation of bound magnetic 
polaronsAi 2 .. The bound magnetic polaron (BMP) is a 
local ordering of magnetic moments induced by the ex- 
change interaction with a localized carrier. Most BMP 
studies, both experimental and theoretical, have been 
performed so far in a regime of the mean-field approach 
(when N Mn a 3 B > 1 Here, N Mn = xN is the 

concentration of magnetic ions, Nq is the concentration 
of cations, and cib is the localization radius. This ap- 
proach implies that the localized carrier interacts with a 
infinite number of Mn 2+ ions, and hence it is typically 
applicable for high and moderate manganese concentra- 
tions (x > 0.01). At present systems with a countable 
number of magnetic moments attract growing interest^. 
This corresponds to the very dilute regime NMn a % ^ 1; 
i.e. when on average one carrier interacts with one (or 
few) Mn ion(s). Clearly, at this condition the mean-field 
approach is not an appropriate model. 

Up to now bound magnetic polarons have not been 
observed in a very dilute regime. This is because their 
observation requires the polaron formation time tf to 
be shorter than the life time r of nonequilibrium carri- 
ers, Tp < t. In DMS with x > 0.1 this condition is 
fulfilled. However, in samples with lower Mn concen- 
trations x < 0.01 the polaron formation time is of the 
order t f ~ 10~ 4 10~ 6 si2.il. Therefore, the BMP for- 
mation process is interrupted by exciton recombination 
(with characteristic time r ~ 10 -9 s)i£. 

In this paper we study BMP in bulk Cdi_ x Mn x Te 
samples with very low Mn concentration (x < 0.01), 
i.e. at the condition Nu n a% « 1. In order to avoid 
the above problem of BMP detection in the very dilute 
regime we use resonant (selective) excitation of donor- 
acceptor (DA) pairs with transition energies^ 



Here, Ed and Ea are the binding energies of donors and 
acceptors respectively, E g is the band gap and e is the di- 
electric constant in Cdi_ x Mn x Te. As follows from Eq. 
by changing the excitation energy E ex = Ed a one can 
selectively excite pairs with a given DA separation p. Be- 
cause of the weak overlap [~ exp(— p/clb)} of the donor 
and acceptor wavef unctions, the pairs with large p are 
characterized by a long recombination time r > t f . As 
a result the bound magnetic polaron forms before the re- 
combination occurs. Without polaronic effects such reso- 
nantly excited pairs should manifest themselves in emis- 
sion spectra as a narrow line of width within either the 
inverse recombination time (H/t) or the laser line width. 
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The bulk Cdi_ x Mn x Te samples with Mn concentra- 
tions x = 0.003, 0.005, 0.02 have been grown by the 
Bridgman technique. The crystals are nominally un- 
doped. The concentration of residual impurities (both 
donors and acceptors) is about 10 16 cm -3 , which is much 
lower than the Mn concentration Ny[ n . As grown crystals 
are usually p-type due to excess of Te or Cd-vacancies. 
Note, for acceptors with the Bohr radius as = 10 A and 
x = 0.005 (N U n = 7 x 10 19 cm" 3 ) the relationship is 
JVM„a| = 0.07. In the experiments we used 4x4x0.3 mm 
pieces cleaved out along the (110) plane from the mas- 
sive monocrystals. All measurements were carried out 
at a temperature T — 1.6 K. External magnetic fields 
were applied either perpendicular to the sample plane 
(Faraday geometry) or in the sample plane (Voigt ge- 
ometry). The photoluminesccnce (PL) was excited by a 
He-Ne laser or a tunable Ti-sapphire laser pumped by an 
Ar-ion laser. The excitation density was about 5 W/cm 2 . 
The laser beam was directed onto the samples at the an- 
gle close to the normal (axis [110]), and the emission was 
registered in the backscattering geometry. A polarizer (in 
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FIG. 1: (Color online) (a) Zero-field PL spectrum of the 
Cdo.ggsMno.oosTe sample excited by a He-Ne laser (E ex — 
1.959 eV). (b) PL spectra under quasiresonant excitation, the 
excitation energies (E ex ) being shown in the panel. Arrows 
indicate emission bands which are related to the BMP forma- 
tion. PL in the same spectral range but under nonresonant 
excitation [as in (a)] is given in the upper part for comparison, 
(c) PL spectra of DA pairs recorded in a + and a~ circular 
polarizations in an external magnetic field B — 7 T applied 
in Faraday geometry. The excitation conditions are the same 
as in (a), (c) Degree of circular polarization vs magnetic field 
detected at the DA band. 



combination with a quater-wave plate) was (were) used to 
linearly (circularly) polarize the excitation. The degree of 
PL polarization was detected with use of a photo-elastic 
modulator and a two-channel photon counter. 

We found qualitatively very similar behavior in all 
samples studied. In the following we present experi- 
mental results for one piece cleaved out from the bulk 
Cdo.ggsMno.oosTe sample. Figure QJa) shows the PL 
spectrum excited above the CdTe band gap E g (Hc-Nc 
laser with E ex = 1.959 eV). The high-energy line (X) 
is attributed to the bound exciton. At low-energy part 
of the PL spectrum we observe a broad band associated 
with recombination of DA pairs and its phonon replicas. 
These replicas are separated by 20.5 meV, which corre- 
sponds to the LO-phonon energy in CdTei^. Application 
of external magnetic fields in Faraday geometry results in 
in the circularly polarized emission already in low fields 
[Fig. [He)]. Such a behavior is typical for magnetic semi- 
conductors. Figure QJd) demonstrates the degree of cir- 
cular polarization vs magnetic field detected at the DA 
band. It saturates rapidly on a 75% level. 

In order to study bound magnetic polarons we use the 
site selective spectroscopy 14 . We resonantly excite DA 
pairs and record BMP spectra shifted towards lower ener- 
gies with respect to the excitation energy E ex (the Stokes 
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FIG. 2: (Color online) (a) BMP spectrum obtained under 
resonant excitation (Ex.) to the DA band (E ex — 1.551 eV) in 
zero magnetic field, (b) BMP and BMP-polarization spectra 
obtained under the same excitation conditions as in (a) but 
in magnetic field B — 7 T applied in Faraday geometry, (c) 
degree of circular polarization vs magnetic field detected at 
different energies, (d) BMP spectrum obtained in magnetic 
field B = 7 T applied in Voigt geometry when the polarization 
axis of the excitation (E) is parallel to the field direction 
(E||B). 



shift). A set of PL spectra for different E ex is presented 
in Fig.QJb). A new emission band (indicated by arrows), 
being observed on the top of the DA photoluminescence, 
clearly follows the laser line which is a fingerprint of the 
BMP formation. The most pronounced results have been 
obtained for E ex = 1.551 eV, as also shown in Fig.O In 
order to avoid scattered light from the laser we detect 
such BMP spectra at the first phonon replica. 

In zero magnetic field the BMP spectrum represents 
itself as an asymmetric band with the 5-meV-tail towards 
low energies [Fig. [Ha)] . Our estimations (given below) 
show that this value is in agreement with the BMP energy 
as well as with its root- mean-square fluctuations. Such 
a behavior is an essential deviation from the case of high 
Mn concentrations (NmuO.% 3> 1), where the Stokes shift 
is larger than the line width, giving rise to a distinct 
emission line shifted from the excitation by the BMP 

14 

energyii. 

The most drastic changes, however, occur in external 
magnetic fields applied in Faraday geometry [Fig. [3ja)] . 
A new line appears at the resonant energy (i.e., there 
is no Stokes shift at all)^. With growing magnetic field 
the amplitude of this line increases and the width reduces 
greatly. In strong enough magnetic fields (B > 3 T) the 
width saturates on a 240- fj,eV- level [see also Fig. |2Jb)]. 
The circular polarization of this line is about 65%, and 
its dependence on the magnetic field [Fig. [DJc)] follows 
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FIG. 3: (Color online) Evolution of BMP spectra with grow- 
ing magnetic field: (a) experimant and (b) the Monte-Carlo 
simulation. Note, in (b) no fitting parameters are used. The 
PL spectra are shifted for clarity. 



that for the DA emission under nonresonant excitation 
[Fig. Old)]. One more lineament of this new line: Its 
intensity reduces under excitation by the light being lin- 
early polarized along the magnetic field direction applied 
in Voigt geometry [see Fig.^d)]. However, a deviation of 
the field direction from the sample plane or the polariza- 
tion axis from the field direction restores the amplitude 
of this line. 

As we show in Sec. lIII| the evolution of PL spectra pre- 
sented in Fig. [3]^a) cannot be described in terms of the 
mean-field theory, and hence this requires an alternative 
approach. First, we discuss the asymmetric band ob- 
served in zero magnetic field. Resonant excitation gener- 
ates DA pairs with a certain energy Eda = E ex depend- 
ing on donor-acceptor separation p [pairs (i) and (iii) in 
Fig. Ufa)]- DA pairs with transition energies different 
from the excitation energy Eda ^ E exc are not excited 
[pair (ii) in Fig. Ufa)]. The observation of the asymmet- 
ric band in Fig. [2a) we attribute to the BMP forma- 
tion accompanied by relaxation to the ground state. The 
dominant contribution to the polaron shift arises from 
the hole. This is because the absolute value of the p-d 
exchange constant is larger than the absolute value of 
the s-d exchange, and the Bohr radius of the acceptor is 
smaller than that of the donor— In this approach the 
line width is explained by the dispersion of the polaron 
shifts induced by static fluctuations of Mn concentration. 
For instance, pair (iii) in Fig.[4ja) shows a larger polaron 
shift than pair (i). Because on average less than one 
Mn falls into the volume restricted by the acceptor Bohr 
radius, the root-mean-square fluctuations appear to be 
comparable with the polaron shift. Its value can be esti- 
mated from the line width in Fig. [2^a), T ~ 5 meV. 

Strong enough magnetic fields suppress the polaron 
formation 16 , and one would expect the disappearance of 
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FIG. 4: (Color online) Formation of BMP spectra under se- 
lective excitation (Ex.) in (a) zero and (b) strong magnetic 
fields. The roman numerals in brackets enumerate different 
DA pairs. 



the low-energy tail. However, our experiments show the 
opposite behavior: This tail preserves and a new narrow 
intense line appears. We suggest the following expla- 
nation. Acceptor (donor) states are split in an exter- 
nal magnetic field. The photon absorption by the fixed 
DA pair leads to the population of one of the Zeeman 
sublevels [Fig. Hfb)]. Relaxation from the higher-lying 
sublevels to the lowest sublevel [pairs (v) and (vi) in 
Fig. HJb)] results in the low-energy tail. The character- 
istic width of this tail does not change with magnetic 
fields because the Zeeman splittings disperse in the same 
manner as the polaron shifts, i.e. due to the static fluc- 
tuations of the Mn concentration. In case when the low- 
est Zeeman sublevel of DA pairs is excited [pair (iv) in 
Fig. QJb)] a further relaxation is not possible, and the 
narrow line is observed at the excitation energy. With 
growing magnetic fields the width of this line should de- 
crease and saturate in fields when the magnetization of 
magnetic ions saturates (B > 3 T). 

It is now clear why the narrow peak is suppressed when 
the polarization axis of the excitation light is exact par- 
allel to the magnetic field direction: Due to optical se- 
lection rules the excitation to the ground state [pair (iv) 
in Fig. Hfb)] is forbidden in this configuration. How- 
ever, even a small deviation from this geometry restores 
the perpendicular-to-field component of the polarization 
vector, and the peak rises up. 



III. THEORY 

In order to calculate the PL spectra under selective 
excitation one needs to find the distribution function of 
the Stokes shift induced by the polaron formation in zero 
magnetic field. On the other hand, in case of strong 
magnetic fields the distribution function is determined 
by the Zeeman splitting dispersion. The hole (electron) 
spins interact with Mn ions via the p-d (s-d) exchange 
interaction. In case of localized carriers this interaction 
can be written as 

Hsd = aJ2~& -Sl^Kn), (2) 

i 
i 

where ~sl and J are the electron and hole spin opera- 
tors, 4" e an d are their wavefunctions, respectively. Si 
is the spin operator of the i-th Mn ion, and is a radius- 
vector assigning Mn location in the crystal. The ex- 
change constants in Cdi_ x Mn x Te are a = 14850 meVA 3 
and (3 = 59400 meVA 3 . The donor and acceptor Bohr 
radii are equal approximately 60 A and 10 A, respec- 
tively. This implies that the spin splitting of the Mn ion 
removed by the Bohr radius from the acceptor (donor) is 
1.28 meV (1.5 /ieV). In other words, at a temperature 
T = 1.6 K (k B T = 0.14 meV) the electron contribu- 
tion can be neglected, and all Mn spins in the vicinity 
of the acceptor are aligned along the exchange field of 
the hole. The hole exchange energy in the ground state 
is Eh = j/3Y^i^h( r i)- Its mean value averaged over 
all Mn configurations is (Eh) — jflNoX, where Nq = 
0.015 A -3 is a concentration of cations. For x — 0.5% it is 
equal (Eh) = 5.0 meV. The root-mean-square fluctuation 
SE h = y/(El) - (E h ) 2 is equal SE h = 5.9 meV. And the 
mean value of the hole energy in the fluctuating fields of 

nonpolarized Mn ions SE h f = \f3\J 'n x^^- J l^rfV 
is equal SEhf ~ 5.5 meV. It follows from these estima- 
tions that the fluctuations of the exchange energy and its 
mean value are nearly the same. Therefore the Gaussian 
statistics and, as consequence, the mean-field approach 
is not applicable. 

In case of low Mn concentrations it is difficult to ob- 
tain an analytical solution for the Stokes shift distribu- 
tion function. Therefore, in order to calculate the PL 
spectra under selective excitation we use numerical com- 
puting based on the Monte-Carlo method. First, a ran- 
dom distribution of magnetic ions in a sphere with a ra- 
dius equal to ten acceptor Bohr radii lOa^ is generated 
[see also Fig. [5jb)] 17 . Each Mn ion has a random spin 
orientation in accordance with the equilibrium distribu- 
tion for a given temperature and magnetic field strength. 
Then, the total exchange field acting on the hole bound 
to the acceptor placed in the center of the sphere is calcu- 
lated. After that, the eigenstates of the hole in this total 
exchange field and the transition probabilities on these 




FIG. 5: (Color online) (a) Calculations of the PL spectra by 
the Monte-Carlo method and within the approach considering 
the nearest-Mn ion. Note, the nearest-Mn approach gives the 
same result for zero and strong magnetic fields, (b) Random 
distribution of acceptors (A), donors (D) and Mn ions. 

states are evaluated. Polaron energies are determined for 
each generated Mn distribution by adding the exchange 
filed of the hole acting on Mn ions to the external mag- 
netic field. The Stokes shift is calculated as an energy 
difference between the initial and final states. Finally 
the distribution function of the Stokes shift is obtained 
after summation over 50000 realizations. 

The calculated PL spectra for different magnetic fields 
are shown in Fig. 02L). Not only qualitative but also 
quantitative agreement with experimental data [compare 
with Fig. Ufa)] is achieved. Note, no fitting parameters 
are used in the calculation. 

Because of low Mn-ion count in the vicinity of the ac- 
ceptor (on average there are 0.3 Mn ions inside the sphere 
of radius as) it is natural to assume that the nearest 
Mn ions give dominant contributions to the Stokes shifts 
[Fig. [3b)]. We now consider this approach in details. In 
an external magnetic field the hole states are split in four 
sublevels characterized by the different projection J z of 
the hole angular momentum on the field direction. One of 
them with J z = 3/2 is a ground state. The contribution 
to the PL spectra from this state is a delta function 6(e) 
(if neglecting native homogeneous broadening) as shown 
in Fig. BJb) [pair (iv)]. Excitation in one of three other 
states is followed by relaxation to the ground state [pairs 
(v) and (vi) in Fig. HJb)]. Therefore, the shift of the 
emission energy with respect to the excitation energy is 
described by the distribution function of the nearest (to 
the hole) Mn ion, Fjvjv(e)- Then, under unpolarized ex- 
citation the distribution function F(e) of the Skokes shift 
e is written as 

F(e) = W 3/2 5(e) + W 1/2 F NN (3e/2e Q ) + 

+W 1/2 F NN (3e/4e ) + W 3/2 F NN (e/2e ) , (4) 

where W 3 / 2 = 3/4 and Wi/ 2 = 1/4 are the probabili- 
ties to excite the hole with J z = ±3/2 and J z = ±1/2, 
respectively. And the characteristic exchange energy is 
£o = |/3tf£(0) = 23.6 meV. The probability dP NN (r) 
that the nearest Mn ion is located in the spherical layer 
(r,r + dr) [see Fig. [SJb)] is given by the probability to 
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find one Mn ion in this spherical layer being equal to 
xN Q A-Kr 2 dr, and the probability that there are no Mn ions 
inside the sphere of radius r being equal to Pq (r) . Here, 
Po(r) is the Poisson distribution P n (n) = exp(— n) 
with n = and n = ^Trr 3 xNo. Then one writes 

dP NN (r) = 4nr 2 xN exp (-47rr 3 xN /3) dr . (5) 

Recall, the energy of the exchange interaction between 
the hole and Mn ion depends on a distance r as ^\(r). 
i.e. e(r) = Eq exp(— 2r/as) and dr = —aB/(2e)de. Fi- 
nally, one obtains the distribution function vs energy e 
as -F/v n = APn n /de 
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(6) 



with 7 = (ir /6)xN a 3 3 . 

Figure EJa) shows that the nearest-Mn approach de- 
scribes well the PL spectra in high magnetic fields. How- 
ever, the agreement becomes worse for small e in zero 
magnetic field, while the high-energy tail is still described 
quite well. The reason is that for the complex 'hole + 
nearest-Mn ion' the relaxation is not possible at the exci- 
tation to the ground state, and the Stokes shift is absent. 
This results in a (5-function peak appearing in zero as 
well as in high magnetic fields [first term in Eq. Q]. In 
real the contribution to the Stokes shift for e <C £o arises 
from several distant (r 3> ag) Mn ions. In weak mag- 
netic fields the spins of these distant ions are disordered, 
and their orientation induced by the exchange field of 



the hole results in the energy lowering. In strong mag- 
netic fields the spins are already oriented at the moment 
of excitation, further relaxation does not occur, which 
manifests itself as the appearance of the narrow peak at 
e = 0. 



IV. 



CONCLUSIONS 



With use of the site selective spectroscopy of residual 
DA pairs in bulk Cdi_ x Mn x Te samples we are able to ob- 
serve the formation of bound magnetic polarons in a very 
dilute (x < 0.01) regime. Under selective excitation the 
PL spectra show an asymmetric band with low-energy 
tail and the additional very narrow peak inflaming in 
external magnetic fields. This behavior cannot be de- 
scribed in terms of the mean-field approach, frequently 
applied to diluted magnetic semiconductors (DMS). Our 
findings suggest that the low-energy tail is formed due 
to the exchange interaction of the hole (bound to the ac- 
ceptor) with the nearest Mn ion, and is described by the 
Poisson distribution. Unexpectedly, the dominant con- 
tribution to the narrow peak arises from the interaction 
with several Mn ions remote from the hole. 
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